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O ■ ABSTRACT 
u: We report the results of a pilot study with the Expanded Very Large Array (EVLA) of 12 CO 

fH ' J = 1 — emission from four submillimetre-selected galaxies (SMGs) at z = 2.2-2.5, each 

O | with an existing detection of 12 CO J — 3 — 2, one of which comprises two distinct spatial 

components. Using the EVLA's most compact configuration we detect strong, broad (medi- 
ans: 990 kms -1 FWZI; 540 km s -1 FWHM) J = 1 — line emission from all of our targets 
- coincident in position and velocity with their J = 3 — 2 emission. The median line width 
ratio, <ti_o/<73-2 = 1.15 ± 0.06, suggests that the J — 1 — is more spatially extended than 
the J — 3 — 2 emission, a situation confirmed by our maps which reveal velocity structure in 
several cases and typical sizes of ^16kpc FWHM. The median brightness temperature (lb) 
C*~) ■ ratio is r 3 _ 2 /i-o = 0.55 ± 0.05, consistent with local galaxies with Ljr > 10 11 L Q , noting 

that our value may be biased high because of the J = 3 — 2-based sample selection. Naively, 
this suggests gas masses roughly 2x higher than estimates made using higher- J transitions 
of CO, with the discrepency due entirely to the difference in assumed Tb ratio. We also esti- 
mate molecular gas masses using the 12 CO J = 1 — line and the observed global Tb ratios, 
assuming standard underlying Tb ratios for the non-star-forming and star-forming gas phases 
as well as a limiting star-formation efficiency (SFE) for the latter in all systems, i.e. without 
calling upon Xco (= £*)• Using this new method, we find a median molecular gas mass of 
(2.5 ± 0.8) x 10 10 M , with a plausible range stretching up to 3 x higher. Even larger masses 
cannot be ruled out, but are not favoured by dynamical constraints: the median dynamical 
mass within R ~ 7kpc for our sample is (2.3 ± 1.4) x 10 11 M Q , or ~6x more massive 
than UV-selected galaxies at this epoch. We examine the Schmidt-Kennicutt (S-K) relation 
for all the distant galaxy populations for which CO ,7=1 — 0orJ = 2 — 1 data are avail- 
ed ■ able, finding small systematic differences between galaxy populations. These have previously 

been interpreted as evidence for different modes of star formation, but we argue that these 
differences are to be expected, given the still considerable uncertainties, certainly when con- 
sidering the probable excitation biases due to the molecular lines used, and the possibility of 
sustained S-K offsets during the evolution of individual gas-rich systems. Finally, we discuss 
the morass of degeneracies surrounding molecular gas mass estimates, the possibilities for 
breaking them, and the future prospects for imaging and studying cold, quiescent molecular 
gas at high redshifts. 

Key words: galaxies: evolution — galaxies: high-redshift — galaxies: starburst — infrared: 
galaxies — radio lines: galaxies 
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1 INTRODUCTION 

The star-formation density contributed by ultraluminous infrared 
galaxies (ULIRGs) appears to increase out to z £ 2 (e.g. Chap- 
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man et al. 2005; Wardlow et al. 2010). These galaxies are often 
heavily obscured by dust and hence the bulk of their luminosity is 
radiated in the rest-frame far-infrared (far-IR) and observed in the 
submillimetre (submm) band, hence their epitette: 'submm galax- 
ies' (SMGs). SMGs have the potential to form the stellar mass of 
an L* galaxy in a single event (e.g. Lilly et al. 1999; Smail et al. 
2004; Swinbank et al. 2006; Hainline et al. 2010). To accomplish 
this feat, however, SMGs must have sufficiently large reservoirs of 
cold gas, ^lO 11 Mq . The first searches for molecular gas emission, 
via the J = 3— 2 or 4—3 transitions of 12 CO, were successful in de- 
tecting significant quantities of gas in several SMGs (Frayer et al. 
1998, 1999; Ivison et al. 2001). Subsequent, larger surveys with the 
Plateau de Bure Interferometer (PdBI) greatly expanded this work 
(e.g. Downes & Solomon 2003; Genzel et al. 2003; Kneib et al. 
2005; Greve et al. 2005; Tacconi et al. 2006; Bothwell et al. 2010). 

Observations of the molecular gas within high-redshift galax- 
ies provides powerful insight into the physics of star formation in 
these systems (Solomon & Vanden Bout 2005), and allows com- 
parisons with local systems. More specifically they allow us to: a) 
probe the mass and extent of the reservoir of molecular gas avail- 
able for fueling their prodigious starbursts; b) determine the dy- 
namical mass of the host galaxy, free from the uncertainties arising 
from outflows and patchy dust extinction which plagues optical and 
near-IR spectroscopic studies (e.g. Swinbank et al. 2006; Ivison 
et al. 2010a); c) derive their gas-mass fraction, M gas (H2)/M d yn. 
The dynamical mass is a strong indicator of the 'end' state of such 
systems in the present Universe while a well-determined gas-mass 
fraction indicates their likely evolutionary status at the look-back 
time where they are observed. Finally, gas-consumption timescales, 
(Tgas) = M(H2)/SFR, where SFR is the star-formation rate, 
gives the minimum timescale for the conclusion of their star- 
forming (SF) episodes. 

While the PdBI studies of SMGs represent a considerable ad- 
vance, they are fundamentally limited by their focus on high-J 
12 CO lines whose high excitation requirements (n cr it ~ 10 4 - 
10 5 cm~ 3 , E u /kB ~ 50-150 K) confine the emission from such 
transitions to regions of active star formation, rather than tracing 
the total available reservoir of gas within a galaxy. In the nearby 
archetypal starburst, M 82, such lines would reveal only the highly 
excited, SF molecular gas in its inner 400 pc (Mao et al. 2000; WeiB 
et al. 2001) rather than the more massive, low-excitation gas com- 
ponent that extends ~ 1.7 kpc beyond its centre (Walter et al. 2002; 
WeiB et al. 2005). This also suggests potential spatial biases in the 
high-J 12 CO emission - which exhibit typical half-light radii of 
0.8-2.8 kpc for SMGs (Tacconi et al. 2008) - means they may not 
trace the true dimensions (or kinematics) of the total molecular gas 
distribution. 

A recent comparison of 1 CO line ratios between local, IR- 
luminous galaxies and distant SMGs has provided strong indica- 
tions that several high-redshift systems must contain significant 
amounts of colder, non-SF gas or - for a few compact, extreme star- 
bursts - that they may suffer high optical depths at short submm 
wavelengths due to dust (Papadopoulos et al. 2010). Thus their 
high-J 12 CO line emission may not be a good tracer of the total 
CO-rich molecular gas mass, its distribution, or the total enclosed 
dynamical masses, as is usually assumed (e.g. Tacconi et al. 2006). 
Indeed the low brightness-temperature (Tb) ratios of 12 CO J = 7—6 
(or 6—5) to 3—2 (r 7 _ 6 /3_2 < 0.3) measured in several SMGs (Tac- 
coni et al. 2006) are not typical of dense, warm, SF gas providing 
strong but still circumstantial evidence that high-J 12 CO line stud- 
ies may miss a critical gas component in SMGs. 

The earliest evidence of low, Milky Way-type, global CO line 



excitation in a distant SF system was uncovered in the submm- 
bright extremely red object, HR 10 (Papadopoulos & Ivison 2002). 
Since then a handful of SMGs have been observed in 12 CO J= 1— 
(Greve et al. 2003; Hainline et al. 2006; Ivison et al. 2010a; Swin- 
bank et al. 2010; Carilli et al. 2010; Harris et al. 2010; Frayer 
et al. 2010) and several show evidence of substantial reservoirs 
of 12 CO J = 1 — with r3_2/i-o ~ 0.5 (see also Danielson 
et al. 2010). A systematic under-estimate of molecular gas mass 
via high-J CO line emission could lie at the heart of the appar- 
ent discrepancy between the gas-depletion timescales for SMGs, 
40-100 Myr, based on high-J observations (Greve et al. 2005), 
and the proposed lifetimes of these luminous starbursts, ~300 Myr 
(Swinbank et al. 2006, 2008), though strong feedback events will 
punctuate the evolutionary path of any starburst and these are ex- 
pected to lengthen any putative gas-consumption timescales. This, 
in turn, influences our understanding of the evolutionary links be- 
tween SMGs and other high-redshift populations, e.g. quasars or 
passive galaxies (Chapman et al. 2005; Wardlow et al. 2010) which 
rely on duty-cycle and space-density arguments, which then influ- 
ences our interpretation of their likely descendants at z ~ (Swin- 
bank et al. 2008). Similarly, the spatial bias in the high-J 12 CO 
dynamical tracers employed for SMGs means we might be missing 
important signatures in the velocity fields at large radii, perhaps re- 
flecting disk-like rotation or even evidence for cold-flow accretion 
(e.g. Dekel et al. 2009). 

To address these concerns we have undertaken a pilot study 
with the National Radio Astronomy Observatory's (NRAO's 1 ) Ex- 
panded Very Large Array (EVLA) of 12 CO J = 1 — emis- 
sion from four well-studied SMGs, with Ssso^m ~ lOmJy at 
z — 2.2-2.5, two each in the Great Observatories Origins Sur- 
veys (GOODS) North and European Large Area Infrared Survey 
(ELAIS) N2 submm survey fields (Scott et al. 2002; Borys et al. 
2003), pinpointed accurately via arcsec-resolution radio continuum 
imaging (Ivison et al. 2002; Chapman et al. 2005), with redshifts 
determined via Keck spectroscopy (Chapman et al. 2003, 2005) and 
confirmed via detections of 12 CO J = 3 — 2 at PdBI (Neri et al. 
2003; Greve et al. 2005; Tacconi et al. 2006, 2008). Ultimately, the 
combination of EVLA and PdBI imaging will allows us to compare 
the gas masses, morphologies and dynamics derived from the 12 CO 
J —1 — and 3 — 2 or higher- J lines in these intense starbursts. In 
this first paper we explore those parameters probed by the initial 
phase of our EVLA survey, conducted at relatively low resolution 
using the EVLA's most compact configuration during shared-risk 
time. 

Throughout the paper we use a cosmology with Ho = 
71kms~ 1 Mpc~ 1 , fi m = 0.27, Qa = 0.73 which gives a median 
angular scale of 8.2kpc arcsec -1 for our sample. 



2 OBSERVATIONS AND DATA REDUCTION 

Our sample was chosen such that the 12 CO J = 1 — line is red- 
shifted to the ~33-36 GHz frequency range, where both sub-band 
pairs of the EVLA's new Ka-band receivers can be utilised, yield- 
ing up to 256 MHz of instantaeous bandwidth during the earliest 
shared-risk phase of EVLA commissioning with the new Wide- 
band Interferometric Digital Architecture (WIDAR) correlator (see 
Table 1). We overlapped the two 64 x 2-MHz dual-polarisation 



1 NRAO is operated by Associated Universities Inc., under a cooperative 
agreement with the National Science Foundation. 
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Table 1. SMG sample and observing log. 
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Notes: a CO J = 3 — 2 redshifts from Tacconi et al. (2006); for SMM J123707, we quote the value used to determine velocities in that work. 



sub-bands by a total of 10 channels - to minimise issues with the 
edge channels - and we centred the lines 36 MHz to the red, using 
the redshifts published by Tacconi et al. (2006, 2008). This latter 
precaution proved unnecessary. Our approach ultimately yielded a 
single 236-MHz dataset for each target, ~2,000kms _1 of cover- 
age with 16.7-18.2-km s _1 channels, covering uv ~ 5-100kA. 

Each target was observed for between five and seven 3-hr 
tracks during 2010 April-July (Table 1), a period during which 
~15 functional Ka-band receivers were typically available and the 
EVLA was in its most compact configuration (D), recording data 
every 1 s, with a resultant data rate of ~6.3 GB hr -1 [programme 
AS1013]. The tracks were scheduled flexibly to ensure appropri- 
ate weather conditions for these high-frequency observations. Af- 
ter 2010 May 12, several antennas were plagued by phase jumps, 
but most of these data were salvaged. Three tracks suffered se- 
vere phase jumps or unsuitable weather and were discarded. An- 
tenna pointing was checked every 90 min at 5 GHz. Each track con- 
tained a 5-min scan of 1331+305 (3C286) for absolute flux cali- 
bration (1.87 Jy at 35.1 GHz), and regular (every £5 min) scans of 
the bright, local calibrators, 1302+575 and 1642+394 (~0.35 and 
~7.0 Jy, respectively). The latter has been monitored at 33.75 GHz 
(Davies et al. 2009) and we found consistent flux densities when 
bootstrapping from 0137+331, albeit with evidence of the variabil- 
ity found by Franzen et al. (2009). 

Editing and calibration were accomplished within ALVS 
(31DEC10). For 1331+305, we used an appropriately scaled 22.5- 
GHz model to determine gain solutions; the other calibrators were 
used to determine the spectral variation of the gain solutions (the 
'bandpass'), after first removing atmospheric phase drifts on a 
timescale of 6-12 s via self-calibration with a simple point-source 
model. Despite the significant increase in data volume relative to 
the old VLA correlator, we were able to employ standard AITS 
recipes throughout the data-reduction process. 



3 RESULTS 

3.1 Infrared luminosities 

On several occasions in the subsequent discussion, e.g. §4.3, 4.5 
and 4.6, we shall call upon SFRs determined from rest-frame 8- 
1,000-^im luminosities, Lir, so we start by describing these mea- 
surements. 

For our SMGs, plus 13 others available in the literature with 
measurements of 12 CO J = 1 — or J = 2 — 1 and comparable 
selection biases and general properties (Greve et al. 2003, 2005; 
Hainline et al. 2006; Frayer et al. 2008; Schinnerer et al. 2008; 
Knudsen et al. 2009; Ivison et al. 2010a,b; Harris et al. 2010; Car- 
illi et al. 2010; Swinbank et al. 2010), we derive Lir by fitting 
the spectral energy distribution (SED) template from Pope et al. 



(2008) and that of SMM J2135— 0102 (Swinbank et al. 2010; Ivi- 
son et al. 2010b) to the available photometric data, ignoring the 
radio photometry for those sources known to suffer contamina- 
tion by AGN: SMMJ163650 (Smail et al. 2003), SMMJ02399 
(Ivison et al. 2010a) and SMMJ14009 (Ivison et al. 2000; WeiB 
et al. 2009). The SMMJ2135 SED gave better fits and we adopt 
these values, which are 10-20 per cent higher than those given 
by the Pope et al. (2008) template. We then convert these into 
SFRs, following Kennicutt (1998a). Our five SMGs have a mean 
Lir of (5.4 ± 0.8) x 10 12 L© (see Table 3) and a correspond- 
ing SFR of 930 ± 14OM yr' 1 (for the full sample, Lir = 
(5.5 ± 0.9) x 10 12 L and 950 ± 155 M Q yr" 1 ). 



3.2 Spectra and morphologies 

We construct maps of each source, stepping through the frequency 
range of the data, examining image cubes with velocity resolu- 
tions ranging from 2-16 channels (35-250 km s -1 ). The synthe- 
sised beam for a natural weighting scheme is typically 2.7 by 
2.0 arcsec with a position angle (PA) near —5°, or 3.2 by 2.6 arcsec 
when employing a Gaussian taper that reaches 30 per cent at 80 kA. 
We find strong 12 CO J = 1 — line emission from each of the 
four SMGs, at the expected positions and frequencies derived from 
the 12 CO J = 3-2 emission, though we note that some published 
redshifts/spectra are misleading - e.g. for SMM J123707, Tacconi 
et al. (2006) quote z = 2.490 rather than z = 2.486. We see 
two distinct spatial components in SMM J 123707, so we hereafter 
refer to a sample of five SMGs using the suffixes -NE and -SW 
for SMMJ123707, following Tacconi et al. (2006). We extract a 
spectrum for the emission in an area corresponding to around three 
synthesised beams, centred on each of our targets (SMM J 123707- 
NW and -SW are combined in this case) and show these in Fig. 1, 
after Hanning smoothing with a four-channel triangular kernel. All 
sources display broad lines, with full widths at zero intensity of 
~ 1 ,000 km s ~~ 1 . We measure the line fluxes by integrating the emis- 
sion across the velocity ranges reported in Table 1 (no continuum 
correction is applied as we detect no significant continuum emis- 
sion from any of the SMGs, 3er < 60 /iJy, as expected). We assess 
the uncertainty in the velocity-integrated fluxes from the variance 
of the off-line emission and we report the flux and associated error 
in Table 2. We also list the CO J = 3 — 2 line intensities (from Tac- 
coni et al. 2008) for our five SMGs, having checked for consistency 
with the low-resolution fluxes of Greve et al. (2005). 

We calculate the intensity-weighted redshift for the emission 
in the velocity window and report this in Table 2. It is clear from 
Tacconi et al. (2008) that some of the 12 CO J = 3-2 lines are 
better fit by two Gaussian profiles and this is also true for several 
of the 12 CO J= 1— lines; however, for the sake of uniformity we 
determine the intensity-weighted second moment of the velocity 
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Figure 1. CO J = 1 — spectra and lightly CLEANed images for the SMGs in our study. On the left we show the spectra, integrated over a region of ~3 beam 
areas in the datacubes, with the corresponding PdBI CO J = 3 — 2 spectra from Tacconi et al. (2006) shown as dot-dashed lines, scaled by 9 — 1 X to be on the 
same Rayleigh- Jeans T b scale. Arrows indicate the CO J = 3— 2 line centres. The CO J = 1— spectrum of the NE component of SMM J123707 is shown in 
red. All our targets are well-detected with CO J = 1 — FWHM of ~300-800km s _1 , marginally broader than the corresponding CO J = 3— 2 emission. The 
frequency range used to create the images in the right-hand panels are indicated below each spectrum. A typical ±l-cr error bar is shown at the top right of 
each panel. Velocities are relative to the CO J = 1 — redshifts listed in Table 2 and the spectra have been Hanning smoothed with an four-channel triangular 
kernel with 4-MHz (~35-kms _1 ) spectral resolution. Right: CO J = 1 — emission integrated over the channels indicated in the corresponding spectrum, 
displayed as contours (—3, 3, 4... Xtr where a is typically ~20 /iJy beam^ 1 ), superimposed on linear greyscales of the 1.4-GHz continuum emission. We see 
that all targets are marginally or well resolved at the resolution of our map, with SMM J 123707 comprising two components. The beam is shown in the bottom 
left corner of each image. Insets show 6- by 6-arcsec postage stamps of CO J = 1 — for SMMJ123549 (a 220-km s^-wide bin centred at — 150 kms - 1 , 
displayed as contours, superimposed on a grey scale of the total CO J = 1 — emission) and SMM J 1 23707 (contiguous 390-km s ~ 1 -wide blue and red bins 
represented by the contours and the underlying greyscale, respectively, illustrating the velocity gradient across SMM J123707-NE). 
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distribution and report this as the equivalent 1-u width of a Gaus- 
sian with the same second moment in Table 2. We find a median 
FWHM of 540 ± 110 km s -1 for our target SMGs, which is nearly 
double the median FWHM of 300 ± 20kms _1 measured for the 
12 CO J = 1-0 emission in ULIRGs by Solomon et al. (1997). 
For comparison, we also calculate the equivalent 1-u width of the 
12 CO J = 3 — 2 emission in the same velocity range from the spec- 
tra shown in Tacconi et al. (2008), finding a median ratio between 
the two line widths of cti_o/(7 3 _ 2 = 1.15 ± 0.06. Assuming that 
the dynamics of the central regions of these systems are dominated 
by baryonic gas disks then this suggests that the 12 CO J = 1 — 
is more spatially extended than the 12 CO J = 3 — 2 emission, a 
situation confirmed by our maps. 

Fig. 1 also shows the maps of the 12 CO J = 1 — emission. 
We mark the centroids of the 12 CO J — 3— 2 emission for all of the 
sources, confirming that the 12 CO J = 1—0 emission is spatially co- 
incident in all cases, although we note that in SMM J163658 there 
is evidence that the 1.4-GHz continuum is offset to the north of 
the J = 3 — 2 emission which is, in turn, offset to the north of 
the J = 1 — centroid - a situation reminiscent of the complex 
AGN/starburst, SMM J02399— 1 36 (Ivison et al. 2010a). 

Even though these 12 CO J = 1 — maps were taken with the 
most compact EVLA configuration, they resolve all of the systems 
except SMMJ123707-SW, one of two 12 CO J = 1-0 emitters 
seperated by ~3 arcsec (~23 kpc) in the field of SMM J123707, as 
previously identified in 12 CO J = 3-2 with PdBI and at 1.4 GHz 
with the VLA by Tacconi et al. (2006). The north-eastern compo- 
nent within SMM J 123707 and the other three SMGs have typical 
FWHM of ~2 arcsec or ~ 16 kpc. These sizes are larger than the 
0.8-2.8 kpc half-light radii deduced from the higher- J transitions 
by Tacconi et al. (2008), as expected from the larger line widths. 
Combining our data with the forthcoming C-configuration EVLA 
observations of 12 CO J= 1 — will ensure we have sufficient sen- 
sitivity and resolution to confirm this suggestion. 

We also find evidence for velocity structure within several of 
the sources. In SMMJ123549 we identify a northern spur to the 
J = 1-0 emission (see inset, Fig. 1), centred at — 150 km s _1 , which 
is not seen in J = 3— 2. This feature is coincident with a prominent, 
narrow feature in the J = 1—0 spectrum and suggests that the system 
may contain two components (cf. 4C 60.07 - Ivison et al. 2008). 
We find a velocity difference of just 80 ± 60 km s _1 between the 
two components of SMM J 123707, suggesting that they are orbit- 
ing one another in the plane of the sky with a physical separation 
of ~25kpc. We also see velocity structure within SMMJ123549 
and most strikingly a velocity gradient of ~600kms _1 across 
SMM J123707-NE over a spatial scale of ~2 arcsec (~ 16 kpc - see 
inset, Fig. 1), though our spatial resolution is currently insufficient 
to determine whether the velocity structure is well ordered. 



3.3 Brightness temperature ratios: average gas excitation 

Combining our new 12 CO J = 1-0 data with the 12 CO J = 3- 
2 luminosities from Tacconi et al. (2008), we derive a weighted- 
mean Tb ratio of r , a _ 2 /i-o — 0.56 ± 0.05 (the median ratio and 
bootstrapped error are 0.55 ± 0.05). Including similar observations 
of the same transitions in the five z ~ 2 SMGs from Ivison et al. 

2 The ELAIS N2 continuum image was recreated from the 
raw data [programmes AI91 and AD432] to avoid the astro- 
metric issue described by Morrison et al. (2010). See also 

http://science.nrao.edU/evla/archive/issues/#0 10. 



(2010a), Harris et al. (2010) and Swinbank et al. (2010), we derive 
a weighted-mean Tb ratio of r 3 _ 2 /i_o = 0.65±0.02 and a median 
of 0.58 ± 0.05. The distributions of r 3 _ 2 /i-o for our sample and 
the literature sources are statistically indistinguishable. 

The weighted-mean T b ratio for the Yao et al. (2003) sub- 
sample of local IR-luminous galaxies with rest-frame 8-1,000-pim 
luminosities, Lm, ^ 10 11 Lq, is J" 3 _ 2 /i-o = 0.57 ± 0.06 with a 
median of 0.63 ± 0.10, similar to that of high-redshift SMGs (see 
also Harris et al. 2010). This is also consistent with the average ex- 
citation conditions of local LIRGs studied by Yao et al. but Lm and 
the molecular gas mass is typically an order of magnitude larger for 
the SMGs. This similarity, while reassuring, does not mean that the 
physical state of the molecular gas in SMGs and LIRGs is identi- 
cal. Indeed, r 3 _ 2 /i-o is often low for the LIRGs studied by Yao 
et al., which is atypical of dense, SF molecular gas (where this ra- 
tio would approach unity); this is thought to be caused by a diffuse, 
unbound, warm gas phase that is frequently found in galactic nu- 
clei. The physical scale sampled by our beam in SMGs at z ~ 2.4 
(~20 kpc) is larger than that employed for the LIRGs studied by 
Yao et al. (~2-6 kpc) and may sample an extended, cold molecular 
gas phase that can suppress the global r 3 _ 2/1 _ (an issue discussed 
in detail in §4.3). 

We highlight two features of our r 3 _ 2 /i_o distribution. 
Firstly, we are using what is effectively an 12 CO J = 3— 2-selected 
sample for this analysis and it is possible that SMGs with lower 
r 3 _ 2 /i-o ratios remained undetected by the 12 CO J = 3 — 2 sur- 
veys and were thereby excluded from our study. We note that a pro- 
gramme to search for 12 CO J= 1— emission using the Zpectrom- 
eter instrument on the Green Bank Telescope (GBT) from SMGs 
which were undetected in 12 CO J = 3 — 2 by PdBI has been ter- 
minated due to problems with GBT's Ka receiver. Secondly, we 
stress that there is a wide range in the apparent r 3 _ 2 /i_o values 
within the SMG sample but that formally these are consistent with 
our error-weighted mean value, given the quoted uncertainties, so 
we have yet to uncover the variation in r 3 _ 2 /i_ seen locally (Yao 
et al. 2003). 



3.4 The plausible range of molecular gas masses 

The key advantage of observing the lowest 12 CO J =1—0 transition 
is that it allows us to determine their molecular gas mass and its 
distribution in a manner identical to that used in the local Universe, 
which permits a direct comparison of these systems to local IR- 
luminous galaxies. A firm lower limit on the molecular gas mass 
can be obtained by assuming local thermodynamic equilibrium and 
an optically thin transition where: 

Jg52Uxggr~ 0.08 [f c -^(£PkWp i ))]- 1 

/ [CO/HaA ' 1 m b 

X ^ W-4 ) Kkms-l pc^ ' ^> 

with To = E u /k B ~ 5.5K, J(T) = T (e T ° /T - l)" 1 , 
Tb g = (1 + z) Tcmb = 9.52 k (the temperature of the cos- 
mic microwave background at z — 2.5), gi — 3 (the degener- 
acy of level n = 1), Z ~ 2Tk/T (the partition function), and 
where [CO/H 2 ] = 10~ 4 is the CO abundance in typical molecu- 
lar clouds (or a Solar-metallicity environment - Bryant & Scoville 
1996). Note that Xco is sometimes also known as a and that we 
shall ignore its cumbersome units hereafter. The line luminosity, 
L' coi _ , is the velocity- and area-integrated brightness tempera- 
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Table 2. SMG observed properties. 



IN dlllc 


z CO(l-(>) 


J l-0 


^3—2 


r 3-2/l-0 


r 6-5/3-2 


°l-0 




It 






(Jykms- 1 ) 


(Jy km s -1 ) 






(kms- 1 ) 


(kms- 1 ) 


(kpc) 


SMMJ123549 


2.2015 ± 0.0002 


0.32 ± 0.04 


1.6 ±0.2 


0.56 ± 0.10 


0.36 ± 0.08 


230 ± 20 


1050 


7± 2 


SMMJ123707-NE 


2.4870 ± 0.0005 


0.09 ± 0.02 


0.32 ±0.09 


0.40 ± 0.14 




200 ± 35 


830 


7±2 


SMMJ123707-SW 


2.4861 ± 0.0004 


0.12 ± 0.02 


0.59 ±0.09 


0.55 ± 0.12 




140 ± 30 


920 


< 3 


SMMJ163650 


2.3847 ± 0.0004 


0.34 ± 0.04 


2.3 ±0.3 


0.75 ± 0.13 


0.23 ±0.05 


330 ± 35 


1340 


6±2 


SMMJ163658 


2.4494 ± 0.0002 


0.37 ±0.07 


1.8 ±0.2 


0.54 ± 0.12 


0.33 ±0.07 


295 ± 10 


990 


11 ±3 



Note: a FWHM = 2V21n2 X a. 



ture in the source reference frame, Z/ C01 _ = f Av f A TbdAdV, 
with units Kkms -1 pc 2 (e.g. Solomon et al. 1997): 



Lr, 



3.25 x 10 



I (BSE.)* f [Scoi-odV], (2) 

V ^CO / Jy 



where Di um is the luminosity distance in Mpc, vco is the rest- 
frame frequency (in GHz) of the 12 CO J = 1 — transition and the 
velocity-integrated flux density is in units of Jykms -1 . 

Making no assumption about the state of the molecular gas, 
we set the minimum temperature range to Tki n ~ 15 — 20 K, which 
yields (Xcc> n } ~ 0.45, and this can be used to compute the mini- 
mum plausible molecular gas mass in each SMG, M m i n (for com- 
parison, a temperature range of Tki n ~ 40 — 50 K, typical of SF gas, 
has {Xqq 1 ) ~ 0.65). We determine M max by adopting Xco = 5 
(see §4.1). Note that the contribution from Helium is already in- 
cluded. 

Setting aside issues regarding the appropriate value of Xco, 
it is worth mentioning that the use of CO J = 1 — to estimate 
molecular gas masses in SMGs, rather than J = 3 — 2 or higher- J 
lines, yields ~2x higher masses, simply because we have observed 
that r 3 _2/i-o is significantly below unity. If, for a moment, we 
adopt the so-called 'ULIRG-appropriate' value of Xco, which is 
around 0.8 and has been used widely to estimate molecular gas 
masses in SMGs (e.g. Greve et al. 2005; Tacconi et al. 2006), this 
would yield a median molecular gas mass of (6.1±2.2) x 10 10 M Q 
in these systems. We report the plausible range of molecular gas 
masses in Table 3 and now move on to discuss how best to assess 
the state of the molecular gas, and how we might better determine 
its mass, e.g. locally, the key to choosing an appropriate Xco for 
ULIRGs came from dynamical constraints (Solomon et al. 1997), 
as we discuss in §4.2. 



4 DISCUSSION 

4.1 Physical conditions of the molecular gas 

Two uncertainties plague us when we employ CO transitions to 
trace the mass of molecular gas in high-redshift galaxies, namely: 
1) the assumed Tb ratio, if transitions other than J= 1— are used, 
and 2) the Xco factor, which must be appropriate for the average 
excitation and kinematic state of the molecular gas. The first of 
these has been removed by our observation of the CO J = 1 — 
transition and we now focus our discussion on the second. 

In practice, 12 CO J= 1— emission in giant molecular clouds 
(GMCs) can be optically thick, and M(H.2) = Xco £'coi-o is 
estimated by adopting Xco ~ 5. The latter value of Xco was ob- 
tained from GMC studies in quiescent environments (e.g. Solomon 
et al. 1987; Solomon & Barrett 1991) and is robust (to within 



~2 x ) even in warm and dense SF molecular gas, providing that the 
molecular gas remains reducible to ensembles of self-gravitating 
units in all environments (see Dickman et al. 1986 for the ap- 
propriate formalism). The failure of the latter assumption results 
in Xco ~ 0.3-1.3 in the extreme interstellar medium (ISM) of 
ULIRGs (Downes & Solomon 1998), as the molecular gas reser- 
voir can form a continuous medium rather than an ensemble of viri- 
alised gas clumps. In the literature, Xco = 0.8 is widely adopted 
for SMGs, as noted earlier, but this carries a significant uncertainty, 
even within the ULIRG class, as the original study indicates. Such 
uncertainties can only be exacerbated in SMGs, e.g. more extreme 
velocity fields would act to lower Xco still further, while the pres- 
ence of dominant quantities of dense (n(H2) > 10 4 cm - ' 1 ) SF gas 
could increase Xco - perhaps even back into the Galactic range. 
For these reasons we do not adopt a single Xco value, but explore 
its possible range, aided by the available CO lines, with the J = 1—0 
transition now providing a normalisation to the total molecular gas 
mass. 

Since we want to avoid the assumption that the molecular gas 
reservoirs in SMGs are reducible to ensembles of virialised clumps, 
as well as being able to incorporate any information about the aver- 
age state of the gas (as provided by available 12 CO lines), we adopt 
the following expression: 

Xco = 2.1 ( ^ff ) M (Kkms- 1 pc 2 )" 1 (3) 

where (n(H.2)) and Tcoi-o are the average density and Tt (in 
12 CO J = 1 — 0) of the molecular cloud ensemble (e.g. Bryant & 
Scoville 1996) to be constrained by radiative transfer models of 
available 12 CO line ratios. The expression: 



K v 



(dV/dr) ohs 
(dV/dr). 



1.54 



[CO/H 2 ] 
\/a Aco 



MH 2 )) 
10 3 cm- 3 



-1/2 



(4) 



corrects Xco for non-virial gas motions, where Aco = 
[CO /U 2 ]/(dV/dR) in (kms -1 pc" 1 )" 1 - The parameter, a = 1- 
2.5, depends on the assumed density profile of a typical cloud; here 
we adopt an average value, a = 1.75. Values of K V1I > 1 bring 
Xco towards the lower values reported for ULIRGs (Solomon 
et al. 1997; Downes & Solomon 1998) while K vir < 1 is used 
to exclude dynamically-unattainable kinematic gas states from the 
range of radiative transfer modeling solutions. 

We can then use 1 CO line ratios to constrain the average con- 
ditions of the molecular gas for the SMGs in our sample, and equa- 
tion (3) to obtain the corresponding Xco values. 

For the system with the highest level of global 12 CO excita- 
tion, SMMJ163650, our large- velocity gradient (LVG) modeling 
of r 3 _ 2 /i_o and r 6 _ 5 / 3 _2 yields two plausible solutions, namely: 
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a) T kin ~ 20-30 K with n(H 2 ) ~ 10 4 -10 5 cm~ 3 and K vil ~ 4- 
11, which corresponds to (Xco) ~ 4.5 and b) T kin ~ 90-1 10 k 
with n(H 2 ) ~ 10 3 cm -3 , and K vil ~ 1, with (X C o) ~ 1 for the 
corresponding gas phases. The system with the lowest excitation, 
SMMJ123707-NE, has r 3 _ 2/1 _ ~ 0.4, which is Milky Way- 
like and is also compatible with two ranges of LVG solutions: a) 
T kin = 15 K, n(H 2 ) = 3 x 10 2 cm" 3 , K viT ~ 1 with X C o = 10, 
and b) T kin ~ 25 or 45-100 K, n(H 2 ) = (1-3) x 10 2 cm -3 , 
K viT ~ 4-10 (and even up to 30) with (X C o) ~ 0.8. 

Thus Xco ~ 5-10 is certainly compatible with the global 
12 CO line excitation observed in SMGs, as are values as low as 
Xco ~ 0.4-1. This dual range of LVG solutions and correspond- 
ing Xco factors is a general characteristic of such modeling when 
constrained by a small number of line ratios. In our LVG models we 
always find more parameter space associated with the lower rather 
than the higher values of Xco, but this does not make the latter 
less probable, nor does it help break the aforementioned degener- 
acy (although as we see in the following section, dynamical argu- 
ments favour the former solutions). It is important to note that the 
choice between the two ranges of Xco that could apply in SMGs 
is not between quiescent and SF molecular gas: Xco ~ 5-10 can 
be associated both with cold, non-SF gas as well as with the warm, 
dense gas found in SF regions. Both components can be distributed 
widely in IR-luminous spiral galaxies and their similar Xco fac- 
tors actually ensure the robustness of global molecular gas mass 
estimates in such systems when using a common (usually Galactic) 
Xco value (e.g. Young & Scoville 1991). Low Xco values, on 
the other hand, are found in the warm, diffuse, and typically non- 
self-gravitating gas which - for typical spirals - is confined to their 
nuclei (e.g. Regan 2000), and thus a low Xco applies for only a 
small fraction of their total molecular gas mass. In ULIRGs, low 
Xco factors (with a considerable dispersion) apply for all of their 
molecular gas as the aforementioned gas phase is expected to be 
concomitant with the entire gas reservoir (e.g. Aalto et al. 1995) - 
the result of high pressures in the ISM and tidal molecular cloud 
disruption in merger environments. 

The difference between molecular gas phases with high Xco 
(quiescent or SF) and those with low Xco (found mostly in SF 
environments as a second phase dominating the low- J 12 CO emis- 
sion) is mostly due to the optical depth of the 12 CO J = 1 — line. 
This is large for the former (tl-o > 5) and typically small for 
the latter (ti_o ~ 0.5-2) as a result of higher temperatures, lower 
densities and/or velocity gradients, with K V [ T > 5. The proxim- 
ity of some ULIRGs has allowed several studies to break the de- 
generacies (e.g. Lisenfeld et al. 2000; Hinz & Rieke 2006). When 
using low-J 12 CO lines this is possible with high-resolution imag- 
ing that resolves the gas disks in ULIRGs (Downes & Solomon 
1998), while 13 CO and multi-J 12 CO line measurements (e.g. Mao 
et al. 2000; WeiB et al. 2001) that allow the mean optical depth of 
the 12 CO J = 1 — line to be determined can also achieve this - 
see §4.4. Finally, detailed models of CO line emission from turbu- 
lent gas clouds reveal that CO J = 3 — 2 and higher- J lines are 
very sensitive to the presence of dense, self-gravitating cores; these 
lines are significantly under-luminous with respect to CO J= 1 — 
for diffuse, non-self-gravitating, turbulent gas (Ossenkopf 2002). 
The resulting non-linear dependence of high- J ratios such as CO 
J — 7— 6/J = 3 — 2 (available for several SMGs) on the mean den- 
sity of such a diffuse phase (see Figure 9 of Ossenkopf 2002) makes 
them unsuitable for constraining an unobserved CO J = 1 — line 
luminosity, especially when such a phase exists alongside dense, 
self-gravitating, SF gas, as is the case with ULIRGs. 



4.2 Dynamical constraints on Xco and gas mass 

Table 3 lists the stellar masses for our target galaxies, derived us- 
ing Bruzual & Chariot stellar population models 3 (Hainline et al. 
2010). We caution that these mass estimates are systematically 
uncertain due to the potentially complex star-formation histories 
and dust obscuration within SMGs, as discussed by Hainline et al. 
(2010). We also use the measured FWHM velocity of the 12 CO J = 
1—0 lines and the observed semi-major axes to determine dynamical 
masses, correcting for inclination using (sin 2 i) = 2/3 (following 
Tacconi et al. 2008). Using M dyn = 2.la'{_ R/G we derive a me- 
dian dynamical mass within R ~ 7kpc of (2.3 ± 1.4) x 10 11 M©, 
consistent with the masses of SMGs from their resolved dynamics 
(Swinbank et al. 2006) and roughly 6x more massive than UV- 
selected galaxies at this epoch (Erb et al. 2006), calculated in the 
same manner. 

Combining our minimum estimate of the molecular gas mass, 
(M mi „) = (3.4 ± 1.2) x 10 1(, M© and the stellar masses 
from Hainline et al. (2010), (M») = (1.4 ± 0.3) x 10 11 M©, 
we derive a total baryonic mass, (Mbaryon) = (1.8 ± 0.1) x 
10 11 M Q for our five SMGs, comparable with their median dy- 
namical mass. We also determine a median gas mass fraction of 
(Af m i n )/(Mbaryon) = 0.14 ± 0.02, a median gas-to-stars mass 
ratio of (M min )/(Af») = 0.16 ± 0.02 and a total baryonic frac- 
tion, /baryon = (M bar yon) / (M dyn ) = 0.51 ± 0.10. Adopting the 
Milky Way value for Xco would increase the first two of these 
fractions by a factor of ~10x and would result in /baryon ~ 1.5. 
We note that our median baryonic mass is about 50 per cent higher 
than L* at the present day (Cole et al. 2001) suggesting that SMGs 
are indeed the progenitors of massive galaxies. 

The median baryonic mass does not exceed the median dy- 
namical mass for Xco 3. However, we note that the mass-to- 
light ratios used to derive the stellar masses are not consistent with 
the star-formation histories implied by the observed gas mass, the 
current SFRs and the expectation that we are (on average) seeing 
the SMGs mid-way through their burst phase. If we require that 
SMGs are seen half-way through their burst, that the current SFR 
is 50-per-cent efficient and is sustained for the duration of the burst, 
then we can use the observed dynamical masses, 12 CO J = 1 — 
and rest-frame if -band luminosities to derive self-consistent con- 
straints on Xco, the mass of the stellar population prior to the 
burst, and the duration of the burst. We use a constant-SFR model 
from STARBURST99 (Leitherer et al. 1999) to determine the mass- 
to-light ratio in the ii -band of the burst (which is obscured by a 
foreground dust screen with extinction, Av) and a pre-existing, 
unobscured 1-Gyr stellar population whose combined luminosity 
is required to reproduce the observed if -band flux. A disk-like dy- 
namical model for the molecular gas reservoir in the SMGs yields 
Xco < 2, an expected burst lifetime of £l50Myr, a gas fraction 
of <65 per cent, and moderate extinction for the burst, Av £ 7. In 
this scenario Xco = 0.8 is recovered for a model with a gas frac- 
tion of 25 per cent, a pre-existing stellar mass of 1.4 x 10 11 M 
and a burst duration of 50 Myr with Av ~ 20, which has added 
20 per cent to the stellar mass of the galaxy. However, allowing a 
wider range of dynamical models (e.g. a virialised sphere) removes 
this constraint, allowing solutions as high as Xco ~ 5, and so we 
conclude that it is impossible with current information to reliably 



3 For SMM J123707-NE, we scale the mass given by Hainline et al. (2010) 
for SMMJ123707-SW based on the relative 5.8-jUm fluxes of the compo- 
nents. 
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Table 3. SMG derived properties. 





T/ 

■^coi-o 


M a M b 

J ™min> '"max 


m CQ/SF 


^IR 


JVI bcst ran § e 


ivl * 






(lO^Kkms^^c 2 ) 


(10 10 Mq) 




(1O 12 L ) 


(10 10 M Q ) 


(1O 1O M ) 


(10 10 M Q ) 


SMMJ123549 


7.6 ±1.0 


3.4, 38 


1.3 


5.5 ± 1.2 


2.5-7.5 


21 ± 6 


23 ±4 


SMMJ123707-NE 


2.7 ± 0.6 


1.2, 13 


5.0 


4.1 ± 0.8 


4.9-14.7 


8± 2 


18 ±6 


SMMJ123707-SW 


3.5 ±0.6 


1.6, 18 


1.4 


3.1 ± 0.6 


1.5-3.5 


18 ± 3 


^4 


SMMJ163650 


9.3 ± 1.1 


4.2, 47 


0.33 


7.8 ± 1.5 


2.1-6.3 


14 ±4 


41 ±8 


SMMJ163658 


10.6 ±2.0 


4.9, 54 


1.5 


6.4 ± 0.9 


3.2-9.6 


13 ±3 


58 ±4 



Notes: a) adopting Xco = 0.45 - see §3.4 and §4.1; b) adopting Xco = 5 - see §3.4 and §4.1; 
c) assuming SFE = 500 L , and plausibly up to 3x lower- see §4.3; d) from Hainline et al. (2010). 



constrain Xco from dynamical limits on SMG samples. Neverthe- 
less, it is clear that values of Xco ^ 5 are disfavoured. 



Alternatively, we can use r 3 _ 2 /i_rj to determine the relative 
fractions of cold, quiescent (CQ) and warm, SF molecular gas 

racQ/SF = Mcq/Msf from: 



4.3 The molecular gas in SMGs: two phases 

Larger values of Xco have the potential to radically re-shape our 
view of the structure, gas-consumption timescales and evolution- 
ary state of SMGs. From their observed 12 CO line ratios it is obvi- 
ous that they are not dominated by SF, dense and warm gas where 
r 3-2/i-o ~ 1 an d r 6-5/3-2 ~ 0.8-1, as measured for Orion A- 
type clouds and the central SF regions of nearby starbursts such as 
M82 and NGC253 (e.g. Wild et al. 1992; Bradford et al. 2003). 
For SMGs, the average value of r 3 _ 2 /i-o is lower, with some in- 
dividual sources reaching Milky Way-like values, r 3 _ 2 /i-o ~ 0.4 
and r 6 _ 5/3 _2 ~ 0.23-0.35, well below that of SF gas. 

Recalling our discussion in §4.1, such low 12 CO line ratios do 
not necessarily imply the presence of cold, quiescent gas such as 
that typifying quiescent GMCs in the Milky Way. Diffuse, warm 
and highly non-virial gas, concomitant with the SF phase in SMGs, 
can also suppress their global 12 CO ratios. Such a phase is respon- 
sible for the low average r 3 _ 2 /i_o ~ 0.66 in the nuclei of nearby 
IR-luminous galaxies (e.g. Yao et al. 2003; Leech et al. 2010) where 
its presence has been known for some time (Regan 2000). r 3 _ 2 /i-o 
then drops from its intrinsic starburst value, ~ 1 , as the diffuse phase 
contributes extra 12 CO J = 1-0 emission, while the J = 3 — 2 and 
higher-J lines remain dominated by the SF gas. 12 CO line ratios 
normalised to the J = 1 — line luminosity (e.g. r 3 _ 2 /i_ ) would 
then be lower than those intrinsic to SF gas, while ratios involving 
only high-J lines (e.g. r 6 _ 5 / 3 _ 2 ) would be high, being dominated 
by the SF gas phase. This is not the case for any of the SMGs in our 
sample (Table 2). Moreover, LVG modeling finds no average gas 
state that can adequately reproduce both r 3 _ 2 /i_o and r 6 _ 5 / 3 _ 2 
for any of the SMGs in our sample, with fits of only the high-J 
CO ratios yielding r 3 _ 2 /i_o ~ 0.9-1.2, which is much higher 
than observed. Thus, while a broad two-phase differentiation of the 
molecular gas is certainly apparent in SMGs, it is not of the type 
observed in local starburst nuclei and ULIRGs (with its associ- 
ated low Xco)- Significant masses of non-SF molecular gas remain 
the only viable alternative for the low r 3 _ 2 /i_ and r 6 -5/ 3 -2 ob- 
served in SMGs. This is further corroborated by the fact that the 
systems with the largest spatial extent of the 12 CO J = 1 — rel- 
ative to J = 3 — 2 emission show the lowest r 3 _ 2 /i_o ratios (see 
Fig. 1 and Table 2). 

If we assume conservatively that Xco(SF)~ 
Xco (quiescent) ~ 5 for both the SF and quiescent gas in 
SMGs (in practice Xco for the SF gas can be somewhat lower), 
then we arrive at robust upper limits for their total molecular 
gas masses, M max , as described in §3.4 and listed in Table 3. 



r 3-2/l-0 - r 3-2/l-0 
W-CQ/SF - (CQ) . P) 

r 3 _ 2 /i-o - r 3 _ 2 /i_ 

where r^Syi-o ~ ^ an d r 3-2/i-o ~ are set as typical 
intrinsic ratios of these two gas phases. The total molecular gas 
mass, Mbost = Msf(1 + tugq/sf), can then be deduced us- 
ing Lir (§3.1) if we assume that the SF gas phase in all galax- 
ies has the same intrinsic SFE max = Lir/Msf. Its maximum 
value, ~5OOL Mq 1 , is thought to be the result of an Eddington 
limit set by photon pressure on dust in the molecular gas accreted 
by the SF sites (Scoville 2004; Thompson 2009) while an almost 
identical limiting SFE can result from a cosmic-ray-generated Ed- 
dington limit (Socrates et al. 2008). Such 'maximum SFE' values 
have been observed in diverse places: individual molecular clouds 
around OB star clusters in M 5 1 , for the total molecular gas reser- 
voir of Arp 220, as well as for the HCN-bright gas phase of LIRGs. 
Here we adopt the maximum SFE, noting that the true value will 
likely be smaller by a factor of up to 3 x since gas accretion towards 
sites of star formation may not be spherically symmetric; moreover, 
CO emission may also be included from beyond the natal sites of 
the OB stars (where the Eddington limit is set), which explains why 
the HCN-bright gas phase yields higher values of SFE in ULIRGs, 
closer to the maximum value than those associated with the CO- 
bright gas phase (e.g. Gao & Solomon 2004). 

m CQ/SF an d Mbest are listed in Table 3. Reassuringly, the 
maximum globally-averaged (SFE) - estimated using M m i n (via 
X C o n ) f° r the total gas mass and making no discrimination be- 
tween SF and non-SF molecular gas - is around 160-200 L Mq 1 
for most of our SMGs, with only SMMJ163650 approaching the 
maximum value allowed by our aforementioned arguments. 

The median cold/warm gas mass fraction for the SMGs is 1.4 
(Table 3). Indeed, only in one source, SMMJ163650, do we find 
m CQ/SF < 1. strong evidence that a significant fraction of the 
molecular gas in most SMGs must be cold, low-excitation gas, tak- 
ing no part in the starburst. 4 The large m C q/sF values, along with 
evidence for extended 12 CO J = 1 — emission (this study; Ivison 
et al. 2010a; Carilli et al. 2010), suggests that the main gas reser- 
voir in these systems is more widely distributed than the compact, 
maximal starbursts implied by high-J 12 CO imaging (e.g. Tacconi 
et al. 2006). 

We consider that Mbcst provides the most likely range of the 
total molecular gas masses of the SMGs in our sample, albeit with 
a large uncertainty due to the plausible range of maximum SFE. 
Using this new method, we find a median molecular gas mass of 
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(2.5 ± 0.8) x 1O 1O M , with a plausible range stretching up to 
3x higher. If we combine both components in SMMJ123707, this 
becomes (3.2 ± 0.9) x 10 10 M Q . 



4.4 Breaking degeneracies in excitation and Xco via CO 
observations 

During 201 1-12 the EVLA will offer 2 GHz of instantaneous band- 
width, prior to commissioning of the full 8-GHz WIDAR capa- 
bility. This will allow simultaneous 5 imaging of 12 CO and 13 CO 
,7=1 — 0, each having 4x the velocity coverage utilised to date. 

Sensitive 13 CO J =1-0 imaging of SMGs can break the de- 
generacies of the 12 CO spectral-line energy distribution (SLED), 
allowing us to determine Xco by: a) determining the total molec- 
ular gas mass via an independent, typically optically thin line, or 
b) placing diagnostically powerful lower limits on the 12 CO/ 13 CO 
,7 = 1 — line-intensity ratio, Ri-o- For the 12 CO intensities seen 
in our sample, deep EVLA observations can detect 13 CO 7=1-0 
where Ri-o ~ 3-6, typical for cold Galactic GMCs, where 12 CO 
J = 1 — is optically thick (Scoville et al. 1979; Polk et al. 1988). 
On the other hand, the ability to set lower limits of i?i_o > 10 
can decisively break the degeneracy described in §4.3 in favour of 
the diffuse phase with moderate optical depth in 12 CO 7=1 — 
and low Xco- This is because while Ri-o ~ 5-10 is typical of 
Galactic GMCs and quiescent spiral disks (e.g. Sakamoto et al. 
1997; Paglione et al. 2001), Ri- > 10 is found almost exclu- 
sively in starbursts (e.g. Casoli, Dupraz & Combes 1992; Aalto 
et al. 1995) or galactic centres, where a diffuse, warm, non-self- 
gravitating phase dominates the 7=1 — emission. 

For the SMG with the lowest 12 CO excitation in our sam- 
ple, SMM J123707-NE, our LVG models (for an abundance ratio 
[ 12 CO/ 13 CO] = 60) yield i?i_ > 20 for all solutions corre- 
sponding to diffuse, warm, non-self-gravitating gas, while Ri—o ~ 
9-10 is found for the dense, cold, self-gravitating phase. Similarly, 
for the average r 3 _ 2 /i_o of the sample, and Ri-o = 5, all LVG 
solutions correspond to virial or only slightly unbound kinematic 
states (ifvir ~ 1-3) with Xco ~ 3-6, with the best fit found 
for Tk ~ 15 K and n(H2) ~ 10 3 cm~ 3 - conditions typical for 
quiescent GMCs. For r 3 _ 2 /i_o ~ 0.6 and Ri-o ~ 15, all solu- 
tions compatible with the aforementioned Tb ratio have K V [ T > 5 
and a corresponding Xco ~ 0.65-1.2. For the highest 12 CO 
7= 1—0 flux levels observed in our sample, a -~10x weaker 13 CO 
7 = 1—0 line can be detected at > 10cr by the fully upgraded EVLA 
in ^lOOhr. 

Finally, we mention the possibility that there may be a signifi- 
cant optical depth due to dust in compact ULIRGs with very dense 
gas (n(H2) > 10 5 cm -3 ), even at short submm wavelengths. This 
is suspected for nearby ULIRGs such as Arp220, and for some 
SMGs (Papadopoulos et al. 2010), and creates an additional degen- 
eracy for the interpretation of very low (high- J)/(low- J) 12 CO line 
ratios in extreme starbursts. EVLA observations of the high-density 
gas tracer, HCN 7 = 1 - 0, at 88.632 GHz, along with low- J 12 CO 
transitions and high- J (7 = 6—5 and higher) 1 CO transitions with 

4 Any significant AGN contribution to Ljr would lower our estimates of 
the warm SF molecular gas mass and increase the ticq/sf g as mass frac- 
tions reported in Table 3. 

5 The frequency of 13 CO J= 1-0 differs from that of 12 CO 7 = 1-0 by 
5.070/(1 + z) GHz, or ~ 1 .5 GHz for the SMGs targeted here, though note 
that contiguous placement of sub-band pairs is not mandatory when using 
WIDAR. 
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Figure 2. L IR versus L' CQ for those SMGs, local (U)LIRGs and BzK 
galaxies (see legend) with robust CO J = 1 — (or J = 2 — 1) measure- 
ments. Trend lines are fitted to the ULIRGs plus SMGs (solid) and LIRGs 
plus BzK galaxies (dashed), as well as all the samples combined (dotted). 
The resulting lineal' slopes are given in parentheses. Also shown are the 
mean corrections that would have been required in L' co for the SMGs and 
BzK galaxies if their CO luminosities had been derived using CO J = 4— 3 
or 7 = 3 — 2, which would have resulted in an erroneously steep trend. The 
measurements from this paper are circled. 



the Atacama Large Millimetre Array (ALMA), have the potential 
to break this degeneracy (Papadopoulos et al. 2010). 



4.5 Gas depletion timescales and the evolutionary state of 
SMGs 

The gas-consumption timescale, T gas = A7(H2)/SFR, is the time 
needed for star formation to consume a given molecular gas reser- 
voir, in the absence of feedback effects, and can be estimated using 
the plausible range of values for A/sf and A/best - Thus r gas (SF) = 
Msf/SFR yields the shortest feasible duration for the observed 
SMG bursts, while A/ bcs t/SFR = r gas (SF)(l + ra C Q/ S F) is the 
longest time over which SF can continue if the colder material 
can be involved without delays imposed by dynamics or SF-related 
feedback. 

We find a median gas-consumption timescale, (r gas (SF)) = 
12Myr, which lengthens to 28Myr if the colder gas can become 
involved. This is considerably shorter than naive expectations for 
the lifetime of the SMG phase, ~300Myr (e.g. Swinbank et al. 
2006), but these gas-comsumption timescales can be extended if 
the SFE is below 500 Lq Mq 1 , perhaps up to 3 x longer. Moreover, 
we are neglecting the effect of feedback, which is also expected to 
lengthen the total duration of successive SF episodes in SMGs due 
to the need to re-accrete material towards the typically very com- 
pact SF regions (~ 100-500 pc) expected in merger-driven, gas-rich 
starbursts. 
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4.6 The rest-frame Schmidt-Kennicutt relation in the distant 
Universe: a comparison to the local Universe 

The availability of a significant number of 12 CO J = 1 — line 
luminosities in vigorously SF systems at high redshifts presents an 
opportunity to examine the Schmidt-Kennicutt (S-K) relation be- 
tween their gas reservoirs and SFRs, without resorting to the use of 
high- J 12 CO transitions and assumptions about global 12 CO high- 
.11.1 =1 — line ratios (i.e. assumptions about the average state 
of their molecular gas reservoirs). The latter can create artificial 
offsets from the S-K relation for galaxy populations where only 
high- J 12 CO lines have been observed, if an incorrect 12 CO high- 
.11.1 =1—0 ratio is assumed. Such offsets can also be created by the 
application of inappropriate Xco factors in various galaxy classes. 
For our present short discussion we render the S-K relation only as 
a Lir — L' co relation, postponing our investigation of the physi- 
cal relation, SFR — M(H2), where Xco must be considered, for 
a future paper (Greve et al., in preparation). 

In its original form, the S-K relation was established via H I 
and 12 CO J = 1 — measurements of galaxies in the local Uni- 
verse (Schmidt 1959; Kennicutt 1989, 1998b). The Hi + H 2 gas 
surface densities, E(Hl+H2), are related to the star-formation sur- 
face density, E*, by E* oc [E(HI + H 2 )] K , where n ~ 1.4-1.5. 
In a recent, comprehensive analysis of the H I and H2 distributions 
in nearby galaxies, on sub-kpc scales, Bigiel et al. (2009) found 
a power-law with slope N = 1.0 ± 0.2 in regions where the to- 
tal gas content (H I + H2) is dominated by H2. The latter seems to 
be characterised by gas surface densities ^9M Q pc~ 2 , indicating a 
saturation point for H I. 

Assuming that distant LIRGs are mostly H I-poor, as they are 
in the local Universe, with ^20 per cent of the their total gas in H I, 
we can estimate their total mass using the 12 CO .1 = 1 — line. 
Even for the most H I-rich local LIRGs, omission of H I results in 
a gas-mass under-estimate of <;2x, leaving Xco as the dominant 
uncertainty (see §4.1). 

That SMGs follow the S-K relation in its basic observational 
form, Lir — L' co , was first shown by the 12 CO J = 3— 2 and 4-3 
survey of SMGs by Greve et al. (2005). More recently, CO surveys 
of gas-rich disk galaxies at z = 1-2 (Daddi et al. 2008, 2010; Dan- 
nerbauer et al. 2009; Genzel et al. 2010) have shown that they too 
obey a S-K-like relation, albeit apparently offset from the SMGs 
(lower LmJ L' co ). It has been speculated that this offset reflects 
two modes of global star formation in galaxies: SMGs typify those 
galaxies undergoing major mergers, with intense, highly efficient 
bursts of star formation, and BzKs typify galaxies undergoing a 
more leisurely rate of star formation, at lower efficiency. However, 
this apparent offset in the S-K relation is deduced from a compari- 
son of high- J CO observations from SMGs with a mixture of CO 
J = 3 — 2, J = 2 — 1 and 1 — line observations of disk galaxies, 
and the reality of the apparent offset (and its interpretation as dif- 
ferent global star-formation laws for mergers and disk galaxies) is 
therefore hampered by biases and uncertainties in the underlying 
gas excitation, especially in SMGs. 

With an increased sample of SMGs observed in 12 CO J= 1— 
we can make a relatively unbiased comparison between SMGs and 
the z ~ 1-2 disk population. In Fig. 2 we show Lir (§3.1) ver- 
sus L coi _ , populated with local LIRGs and ULIRGs from Pa- 
padopoulos et al. (in preparation) as well as with those SMGs and 
BzK galaxies with reliable measurements of 12 CO J = 1 — or 
J = 2 — 1 (where J = 1 — is not available, we use J = 2 — 1, 
calculating L' coi _ via r2i = 0.75). We see that the SMGs ex- 
tend the S-K relation to higher luminosities, with a low dispersion. 



Fitting linear relations of the form log Lir = a log L C oi-o + P< 
we determine slopes for the various samples, reporting these in the 
labels of Fig. 2. Taking the SMGs alone, we find a = 0.65 ± 0.13 
and fi — 5.7 ± 1.4; local ULIRGs display a shallower slope, 
a = 0.39 ± 0.11 with /3 = 8.4 ± 1.1, though this is highly un- 
certain as they span only a small range in Lir; for LIRGs we see 
a = 0.67 ± 0.05 and /3 = 5.0 ± 0.5, similar to the SMGs. 

Generally, we find slopes that are significantly below unity 
(a ~ 0.5-0.7) for the various samples. This contrasts with the 
steeper slopes (a ~ 1.1-1.5) reported by other studies of low- and 
high-redshift samples (e.g. Greve et al. 2005; Daddi et al. 2010; 
Genzel et al. 2010). These have typically employed high- J 12 CO 
line luminosities for the high-redshift galaxies and L C oi-o f° r me 
local sources, which may artificially steepen the slope. Iono et al. 

(2009) looked at the S-K relation using 12 CO J = 3-2 for both 
low- and high-redshift galaxies, finding a = 1.08 ± 0.03, similar 
to that derived by Gao & Solomon (2004) using HCN J = 1 - 0. 
These two molecular transitions both trace the dense, warm, SF gas 
phase. A common slope of unity in the corresponding S-K relations 
is to be expected if a near-constant SFE max = Lir/Msf underlies 
star formation in all galaxies, with only the dense gas phase avail- 
able as fuel, and an Eddington-type limit (set by photons or cosmic 
rays) setting SFE max (see §4.3). 

Taking the SMGs and ULIRGs together - they both contain 
extreme SF environments, after all - yields a = 0.56 ± 0.05 and 
/3 = 6.6 ± 0.6, whereas taking the LIRGs and BzK galaxies to- 
gether yields a = 0.72±0.05 and,0 = 4.6±0.5. The BzK galaxies 
do not stand out dramatically from the other low- and high-redshift 
samples. There are only three BzK galaxies with 12 CO J = 1—0 de- 
tections (Aravena et al. 2010), but the situation is unchanged when 
we include the BzK galaxies detected in 12 CO J = 2—1 from Daddi 
et al. (2010), adopting r 2 _i/i_o = 0.75 to calculate L GO1 _ . If, 
however, we compare the BzKs with the SMGs, this time using the 
12 CO J = 3 - 2, 4 - 3 and 5 - 4 lines for the latter (the mean re- 
sulting offset in L' co is shown in Fig. 2), then one might argue that 
SMGs and ULIRGs populate a different sequence to BzK and spi- 
ral galaxies, as did Daddi et al. (2010). Thus an apparent displace- 
ment and/or steepening of the S-K relation between various galaxy 
populations may not reflect a true difference in their respective SF 
modes, but rather the strong excitation biases produced by using 
molecular lines with different excitation requirements. Given these 
biases, the relatively small galaxy samples, and their scatter, the 
evidence for different S-K relations between different galaxy pop- 
ulations is weak, as we shall argue in a forthcoming paper (Greve 
et al., in preparation). 

Adding a last cautionary note, Papadopoulos & Pelupessy 

(2010) have shown that while typical present-day galaxies quickly 
settle into S-K-type relations, this may not be the case for gas-rich, 
metal-poor systems in the early Universe. Such galaxies can spend 
sustained periods with their SFR significantly below or above that 
expected from the local S-K relation. If the same gas-rich, metal- 
poor galaxy can deviate strongly from the S-K relation during the 
course of its evolution, then the use of the S-K relation as a tool for 
differentiating between different SF modes must be re-evaluated. 
Indeed, given the dynamic and non-equilibrium ISM in strongly 
evolving, gas-rich systems (where the global mass fractions of the 
various gas phases, e.g. M S F(H2)/M to tai and M(Hl)/M(H 2 ), 
are strongly time-dependent) there could be a simple, underly- 
ing S-K relation of SFR oc Msf(H2), while the different expo- 
nents recovered in various samples are artifacts of strongly evolving 
MsF(H2)/M to tai fractions, further compounded by the choice of 
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observed molecular lines so with different excitation requirements 
for high-redshift systems. 



5 FUTURE PROSPECTS 

5.1 Blind surveys 

Blain, Carilli & Darling (2004) present the number of blind CO 
J = 1 — line detections expected with 4-GHz-wide observing 
bands, based on Blain et al. (2002), updating the work of Carilli 

6 Blain (2002) to reflect the SMG redshift distribution of Chap- 
man et al. (2005). Their 30-34-GHz band is most easily compared 
with our 2 ~ 2.4 SMGs: scaling to our 236-MHz of instanta- 
neous bandwidth, they predict a source density of ~2.4deg~ 2 at 
the flux levels (§2) to which we are sensitive, > 2 x 10 -22 Wm -2 
or ^O.lSJykms^ 1 . 

It is no great surprise, therefore, that we find no robust (>5- 
(j) detections within the four fields of our pilot survey, each with a 
~70-arcsec FWHM primary beam: sensitivities ranging from ~0. 1- 
0.2 Jykms -1 , covering ~10 -3 deg 2 . 

In Blain et al., the galaxies responsible for J = 1 — emission 
at or above the level to which we are sensitive are extremely lu- 
minous - they have higher Lir, than our SMGs. This suggests the 
predictions may be rather pessimistic; of course, we might also ex- 
pect to see an over-density of sources since our fields are centred 
on relatively massive galaxies (e.g. Stevens et al. 2003). 

Based on the predictions of Blain et al., a blind survey utilising 
WIDAR's full 8-GHz bandwidth, would need to cover ~100x our 
current area, to a similar depth, to provide a significant number of 
robust detections and hence a useful test of the predictions. Our sur- 
vey suggests these requirements can be relaxed significantly; it will 
be interesting to see whether the ~2x deeper observations planned 
for our targets using EVLA's C configuration provide evidence to 
support our suggestion. 

5.2 The study of cold, quiescent molecular gas at high 
redshift 

All current molecular-lines studies of distant starbursts, including 
our own, involve samples selected via large rest-frame far-IR lumi- 
nosities and luminous high- J 12 CO lines (J = 3—2 and higher). It is 
thus possible that SMGs with low r'3-2/1-0 ratios remain mostly 
undetected by the 12 CO J = 3 — 2 surveys and were thereby ex- 
cluded from our study, as noted in §3.3. 

The fact that much of the molecular gas in SF systems may 
not be participating in the starbursts has been well-established by 
studies of local LIRGs, where sensitive low- J 12 CO imaging (e.g. 
Weifi et al. 2004) or submm continuum imaging of dust (e.g. Pa- 
padopoulos & Seaquist 1999; Dunne & Eales 2001; Thomas et al. 
2002) could spatially disentangle the cold and extended molecular 
gas and/or its associated dust from the compact starbursts, which 
otherwise dominate their global SEDs and their molecular SLEDs. 

It has only recently become possible to attempt a similar spa- 
tial separation of ISM components in high-redshift systems, via 
tracers that can remain luminous in the cold, low-excitation ISM - 
e.g. 12 CO J = 1 — and submm continuum emission from dust. 
This explains the previous lack of evidence for extended, low- 
excitation molecular gas around SMGs, but also point to several 
obvious routes forward, via EVLA and ALMA. 

Less well known is the prospect of utilising the two fine- 
structure lines of neutral carbon, 3 Pi-> 3 P at 492.160 GHz and 



P2 — > Pi at 809.343 GHz to trace molecular gas and dynami- 
cal mass (Papadopoulos et al. 2004; Papadopoulos & Greve 2004; 
WeiB et al. 2005), exploiting the full concomitance of C and CO 
in molecular clouds (Keene et al. 1996), the simple, three-level 
partition function, the low optical depth and modest excitation re- 
quirements. Moreover, unlike the luminous [C II] fine-structure line 
at 1.9 THz, the neutral carbon lines are not subject to contamina- 
tion by atomic or ionised gas, tracing solely the molecular gas. 
The [Cl] J = 1 — line remains luminous even for UV-shielded, 
cold (Tk ~ 15 K) molecular gas (e.g. Oka et al. 2001) where the 
[C II] line luminosity is negligable because all carbon is neutral and 
T k < AE U ([C Il])/fc B ~ 92 K. Finally, the [C I] lines are accessi- 
ble in the most sensitive bands (3-7, or 84-373 GHz) of ALMA for 
a much wider redshift range (and range of look-back times) than 
[Cll]-2 ~ 0.3-4.9 for J =1-0 and z ~ 1.2-8.6 for J = 2- 1 
versus z > 4 for [C II]. 

The [C I] 3 Pi — > 3 Po line, in particular, remains well-excited 
in quiescent GMCs (P u //cb ~ 24 K, n cr i t ~ 600 cm~ 3 for low 
Tk), while maintaining a favourable LT-correction with respect to 
12 CO J = 1-0 (Sci/Sco ~ 2.75-5.5 for T b ([C i])/T b (CO)~ 
0.15 - 0.30). The 3 P 2 -> 3 Pi line (E u /k B ~ 63 K, n crit ~ 
965 cm~ 3 ) also maintains a A'-correction advantage with respect 
to both 12 CO J = 1 - and J = 7- 6 (the CO line closest in 
rest-frame frequency), while detecting the pair yields an excellent 
thermometer for molecular gas (both lines are available in bands 3— 
7 of ALMA for 2 ~ 1 . 2^1.9) . Thus [C I] J = 1-0, J = 2-1 imaging 
with ALMA and CO J = 1 — imaging with the completed EVLA 
will constitute the most powerful tools for making inventories of the 
global molecular gas and dynamical gas mass in distant galaxies, 
unbiased by the excitation state of the molecular gas and the extent 
of their star-formation regions. 



6 CONCLUSIONS 

We report the results of a pilot study with the EVLA of 12 CO J = 
1 — emission from a small sample of well-studied SMGs at z — 
2.2-2.5, previously detected in 12 CO J = 3-2 using PdBI. 

Using the EVLA's most compact configuration we detect 
strong, broad (~l,000kms _1 FWZI) line emission from all of our 
targets - coincident in position and velocity with their J = 3 — 2 
emission. 

The median line width ratio, CT1-0/03-2 = 1.15 ±0.06, sug- 
gests that the J = 1—0 emission is more spatially extended than the 
J = 3 — 2 emission, a situation confirmed by our maps which re- 
veal velocity structure in several cases and typical sizes of ~ 16 kpc 
FWHM. With the current spatial resolution we are unable to deter- 
mine whether observed gas motions are well ordered, but we find 
no evidence of large-scale flows of cold gas. 

We find a median Tb ratio of r 3 _ 2 /i-o = 0.55 ±0.05, consis- 
tent with local galaxies with Lir, > 10 11 L Q , noting that our value 
may be biased high because of the J = 3— 2-based sample selection. 
Including five systems with similar luminosities from the literature 
we find a median of r 3 _ 2 /i_o = 0.58 ± 0.05 and see no evidence 
for measureable intrinsic scatter within the sample. 

Using the observed 1 CO J = 1 — line emission, naive es- 
timates of the molecular gas masses are around 2x higher than 
previous estimates based on 12 CO J = 3 — 2 with r 3 _ 2 /i-o = 1. 

We also estimate molecular gas masses using the 12 CO J = 
1 — line and the observed global Tb ratios, assuming standard 
underlying T b ratios for the non-SF and SF gas phases as well as 
a common SFE for the latter in all systems, i.e. without calling 



© 0000 RAS, MNRAS 000, 000-000 



12 Ivison et al. 



upon Xco- Using this new method, we find a median molecular gas 
mass of (2.5 ± 0.8) x 10 M®, with a plausible range stretching 
up to 3 x higher. Even higher masses cannot be ruled out, but are 
not favoured by dynamical constraints: the median dynamical mass 
within R ~ 7kpc for our sample, (2.3 ± 1.4) x 10 11 M Q , ~6x 
more massive than UV-selected galaxies at this epoch. 

We find a median gas-consumption timescale, (r gas (SF)) = 
12 Myr, or 28 Myr if the colder gas can become involved. This is 
shorter than naive expectations for the lifetime of the SMG phase, 
but these timescales can be longer if the SFE is below 500 L Q Mq 1 , 
and we neglect the effect of feedback. 

We examine the S-K relation in Lir — L' co for all the dis- 
tant galaxy populations for which CO J = 1 — or J = 2 — 1 data 
are available, finding small systematic differences between popula- 
tions. These have previously been interpreted as evidence for dif- 
ferent modes of star formation, but we argue that these differences 
are to be expected, given the still considerable uncertainties, cer- 
tainly when considering the probable excitation biases due to the 
molecular lines used, and the possibility of sustained S-K offsets 
during the evolution of individual, gas-rich systems. 

We discuss the degeneracies surrounding molecular gas mass 
estimates, the possibilities for breaking them, and the future 
prospects for imaging and studying cold, quiescent molecular gas 
at high redshift. 

We note in ending that if SMGs (and other high-redshift star- 
bursts) are as extended as our observations suggest (up to 20 kpc) 
then even the shortest possible dish spacings of the EVLA are not 
well-matched to their sizes. However, the smaller ALMA dishes 
(especially those in the ALMA Compact Array), if fitted with 
Band-1 receivers, or GBT, are ideal for studying the critical 12 CO 
J = 1 — emission from these galaxies. 
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